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Abstract. The time-dependent non-crossing approximation is used to study the 
transient current in a single electron transistor attached asymmetrically to two leads 
following a sudden change in the energy of the dot level. We show that for asymmetric 
coupling, sharp features in the density of states of the leads can induce oscillations in 
the current through the dot. These oscillations persist to much longer timescales than 
the timescale for charge fluctuations. The amplitude of the oscillations increases as 
the temperature or source-drain bias across the dot is reduced and saturates for values 
below the Kondo temperature. We discuss the microscopic origin of these oscillations 
and comment on the possibility for their experimental detection. 
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1. Introduction 

Quantum effects are likely to play an increasing role in electronic devices as the physical 
size of their components continues to shrink. Quantum dots and qubits are examples of 
devices where quantum effects play a direct role in their function. In a Single Electron 
Transistor (SET), i.e., a quantum dot coupled to two metallic leads, the conductance 
can be drastically enhanced by the Kondo effect which can occur when the quantum 
dot is populated by an odd number of electrons [HElEllllEliniClIH]- The Kondo effect 
is a quantum-coherent many-body state in which a spin singlet state is formed between 
the unpaired localized electron and delocalized electrons at the Fermi energy at low 
temperatures [9]. 

An important issue for the function of any electronic device is how fast it can 
respond to time-dependent perturbations and bias [101 [H]- Several studies of time- 
dependent response of symmetrically-coupled SET have been performed using a variety 
of methods [l2lll3l[Tlll5lll6l[I7l[T8lll9l|20l[2T]. In the case of a sudden switching of 
the dot level, the transient current has been found to exhibit several distinct timescales 
[22l [23l [24l [25] . The fastest timescale is associated with charge relaxation and the 
other much slower timescales are associated with the formation of a Kondo state, 
i.e. spin dynamics. The detailed evolution of the instantaneous currents following 
a sudden change of the dot level has been shown to depend sensitively on external 
parameters such as source-drain bias, external temperature, dot-lead coupling and 
position of the dot level. These above mentioned time-dependent studies have all been 
concerned with quantum dots which are symmetrically coupled to their leads. While 
the effect of asymmetric coupling on the steady-state conductance of an SET has been 
studied both theoretically and experimentally, [26l |23| to our knowledge, the transient 
transport properties of an asymmetrically coupled dot in the Kondo regime have not 
been investigated previously. 

In this paper, we use a recently developed multi-scale many-body transport method 
to study the effect of asymmetric dot-lead coupling on the transient transport in a 
quantum dot [28]. We show that for a quantum dot asymmetrically coupled to two 
leads with sharp features in their Density of States (DOS), the current can display 
sinusoidal modulations for timescales well beyond the fast charge relaxation timescale. 
The frequency of these sinusoidal modulations is given by the energy difference between 
the Kondo resonance and the DOS feature. The amplitude of the oscillations is 
found to increase with decreasing temperature and source-drain bias and saturate for 
temperatures below the Kondo scale. We attribute this phenomenon to an interference 
effect between the Kondo resonance at the Fermi level of the leads and the conduction 
electrons around the DOS feature. The magnitude of the oscillations depends sensitively 
on the structure of the DOS feature of the leads. We show that these oscillations also 
can occur for leads with a smooth DOS but with finite bandwidth. 
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2. Time dependent current in infinite-U Anderson model 

The SET is modeled as a single spin degenerate level of energy Sdot coupled to leads 
through tunnel barriers, 

o" kaa 

I] [V;(£fcc.)cLaCa + h-C.] , (1) 
kaa 

where c\ {ca) and c^.^^. {ckaa) with a=L,R create (annihilate) an electron of spin a in 
the dot level and in the left(L) and right (R) leads respectively. The n^- and Ukaa are 
the corresponding number operators and are the hopping amplitudes for the left and 
right leads. The Coulomb repulsion energy U is assumed to be sufficiently large that 
double occupancy of the dot level is prohibited. In the following, we will use atomic 
units with h = k-Q = e = 1. 

To deal with large Coulomb correlation, we use the auxiliary boson method for 
the Anderson Hamiltonian. Briefly, the ordinary electron operators on the site are 
decomposed into a massless boson operator and a pseudofermion operator as 

= (2) 

with the constraint that the number of massless "slave" bosons plus the number of 
pseudofermions 

g = 6tfe + ^/t/^ (3) 

must be equal to unity. Self-energies for the pseudofermion and massless "slave" boson 
are then projected onto the physically relevant Q = 1 subspace. 

The Green's functions of the dot levels are calculated using the Non Crossing 
Approximation (NCA) by numerical integration on a multi-scale time grid j^. When 
the dot level Sdot lies well below the Fermi level Sp, the spectral function of the dot 
(local density of states) exhibit two features: a broad Fano-like resonance of full-width 

TUe) = 27iY.{\VL{ek)? + \VR{eu)?)5{e - Su) (4) 

k 

around the dot level, and a sharp temperature sensitive resonance at the Fermi level 
(the Kondo peak), characterized by a low energy scale Tk (the Kondo temperature) 

Tk oc —— exp , (5) 



where D is a high energy cutoff equal to half bandwidth of the conduction electrons and 
Ttot corresponds to the value of Ttot{£) at e = Sp. All energy units in this paper will be 
given in terms of Ttot- 

The current flowing through the SET can be calculated as the difference between 
the currents from the left and right leads as 

Iit) = lL{t)-Init), (6) 
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The most general expression for the net current flowing from the left (right) lead is given 
by |29] 

hiR){t) = -2Im{fjt^ J ^e-(*^-*)r^^)(5) 

,^i:,'^r^Hn^ir)^G<{t,t,) + h^n){e)G^{t,t,)]) (7) 

where the coupling functions T^^ji^ 

rLiR){e) = 2npLiR^{e)VLiR){e)Vl^j,){e), (8) 

depend on the DOS of the leads Pl{r){.^) and the hopping matrix elements in Eq.([l]). 
The quantity ^l{r) represents the time- dependence of the single particle energies in 
the left and right leads. In this paper, we consider the case of a small constant bias 
across the leads and no explicit time and energy dependence of the hopping matrix 
elements, therefore VL{R){e) = VL{R-){ef). We further make the assumption that the 
DOS of the leads are the same, i.e. Pl{R){.^) = p(^)- Thus the coupling functions can be 
parameterized as, 

rL(R)(^) = ri(H)P(^) (9) 

where ^^.(r) are constants given by ^^(/j) = 2t:\Vl(^r){£ f)\'^ and they determine the 
asymmetry of the couplings. In terms of these constants, Vtot = (^l + VR)p{ef). 

In Fig. 1 we show the functions p{e) that will be used to model the leads. The 
Lorentzians used in Fig. Ic and Fig. Id to model sharp DOS features have a width 
equal to 0.002 Vtot- The DOS function p{e) has been normalized for all cases such that 
the bands contain the same number of electrons. The bandwidth of the leads is assumed 
to be 2D with the Fermi energy at e = 0. 



a) 



b) 



-D/2 



-D 





-3D/4 



Figure 1. This figure shows the functions p{e) used to model the leads. Panel a) 
corresponds to a rectangular band, b) to a parabolic band, and c) and d) to a parabolic 
band with a Lorentzian feature located at e=-D/2 and e=-3D/4 respectively 

The expression for the current can then be written as 

I{t)= -2f Llm{f dti{G<{t,ti)h{t-ti)+G''{t,ti)fL{t-ti))) 

J —oo 

+ 2f Rlm{f dt,{G<{t,t^)h{t-t^)+G''{t,t^)fR{t-ti))), (10) 

J —oo 
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where 



and 



.^2^''^"^" ' ^^^^ 

i-D rip „ie{t-ti) 

h{t-t,)= ^p{e) \, (12) 
J-D 2n ]_ + e^\- ~i 

i-D fjp „ie{t-ti) 

In these expressions, V represents the source-drain bias. 

The physical Green's functions in Eq. (11 01) in can be expressed in terms of the 
pseudofermion and slave boson Green's functions using a projection approach discussed 
previously [30]. The final expression for the current is 

I{t) = -2VLlm{f dt,{G<,,^,,{t,t,)B''{1:,,t)h{t~t,) 

- iiiGfseudoit,ti)B<it,,t) + G<,,,,it,t^)B^{t^,t))h{t - ti))) 
+ 2f^/m( /* dti(G'<,„,„(t,ti)i?^(ti,t)/i(t - ti) 

J — oo 

- ^{iG^s,u,oit,h)B<{t,,t) + G<,^,,{t,t,)B^{t,,t))fR{t - ti)))(.14) 

This is the main result of this section and will be used in the calculations presented 
below. 

3. Results 

In this section we will analyze the instantaneous current following a sudden change of 
the dot level from a position at ei=5rtot below the Fermi level where, for the present 
finite temperatures, the Kondo effect will be absent to a position at e2=2Ttot closer to 
the Fermi energy where the Kondo effect will be present. For a parabolic DOS function 
(Fig. lb) with D=9rtot, the Kondo temperature in the final state is approximately 
Tk =0.0016rtot. 

We begin our analysis by investigating the effect of asymmetric coupling to the 
leads. We define the asymmetry factor as the the ratio TL/Ttot where Ttot= + ^r- In 
Fig. 2, the instantaneous conductance, G = I(t)/V, is plotted as a function of time after 
the dot level is switched for various asymmetry factors with small bias. The final steady 
state conductances are in perfect agreement with previous theoretical results [31]. The 
transient short timescale associated with Tfot shown in Fig. 2a is due to the formation 
of the broad Fano-like resonance at 62- The transient increase in the instantaneous 
current is due to the charging of the dot. The steady state current is determined by the 
asymmetry factor. For an asymmetry factor of 1 where the steady state current is zero, 
we have verified the accuracy of our numerical approach by showing that the integrated 
current is equal to the change of the charge of the dot level. 
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Figure 2. Black(dotted), blue(dashed) and red(dot dashed) curves in panel a sliow 
the instantaneous conductance vs. time for rectangular DOS in Fig. la immediately 
after the dot level is switched to its final position for asymmetry factors of 0.95,0.9 
and 0.85 respectively. Ttot is fixed with D^gPtot at T=0.0093rfot and V^Tk- The 
beginning of the oscillations is clearly visible in this panel. Black(dottcd), blue(dashed) 
and red (dot dashed) curves in panel b are the continuation of the ones in the first panel 
in the long timescale for the same parameters. 



Fig. 2b shows the instantaneous conductances for larger times on a magnified 
scale. It is clear from this panel that the current exhibits sinusoidal modulations at 
timescales well beyond Tl or Tn. As we reduce the asymmetry factor, the amplitude of 
these sinusoidal modulations starts decreasing and eventually disappears for symmetric 
coupling. The frequency of the conductance oscillations is equal to the bandwidth D of 
the leads. External parameters such as the energy or width of the dot level, asymmetry 
factor, and ambient temperature and source-drain bias only infiuence the amplitude of 
the oscillations. 

In Fig. 3 we show the effect of the the DOS function of the leads p(e) on the 
conductance oscillations. The oscillations do not look like perfectly sinusoidal functions 
because of our use of a finite time step in the numerical solution of the Dyson equations. 
However, the results are fully converged and for a finer time steps we recover almost 
perfectly sinusoidal oscillations. Figure 3 reveals that the DOS can have a pronounced 
effect on the conductance oscillations. The largest oscillations occur for a DOS with a 
peak feature as in Fig. Ic and Id. For a parabolic DOS function (Fig. lb), conductance 
oscillations are still present but not discernible on the scale of the figure. The frequency 
of the current oscillations is equal to the energy difference between the Fermi level of 
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Figure 3. Effect of the DOS function on the instantaneous conductance. Red(dotted), 
dark blue(dot dashed), black(solid) and green(dashed) curves display the instantaneous 
conductance vs. time in the long timescale for cases a, b, c and d respectively in Fig.l 
when the source-drain bias is equal to Y=Tk with fixed Fjot- AH curves are for 
asymmetry factor of 0.9, D=9Ftot at T^O.OOgSFtot. 

the leads and the feature in the DOS function. For the rectangular and parabolic DOS 
function (Fig. la and lb) where the DOS feature is the band cut-off, the period is equal 
to l/D=l/9 and for the DOS functions depicted in Fig. Ic and Id where the DOS 
feature is the narrow Lorentzian peak, the periods are 2/D and 4/3D respectively. 

We have also investigated other shapes of DOS functions. For instance, for 
a triangular DOS, p{e) = 2(1 + e/D), the conductance displays oscillations whose 
frequency is equal to D very similar to what was obtained for the parabolic DOS in 
Fig. lb. For two sharp features at different energies Di and D2 in the DOS we obtain 
conductance modulations of frequencies proportional to {Di + D2) and {Di — D2). 

Figure. 4 shows the effect of temperature and source-drain bias on the conductance 
oscillations. The figure demonstrates that the amplitude of the oscillations decrease with 
temperature and bias across the dot but their frequency remain unchanged. The figure 
shows two important effects. First, it takes longer for the current and thus the damped 
oscillations to decay as the temperature or source-drain bias is reduced. Second, the 
amplitude of the damped oscillations increases with decreasing temperature or source- 
drain bias and but saturates for temperatures and bias below the Kondo temperature, 
which we estimate to be around Tk =0.0016rtot in these systems. 
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Figure 4. Panel a shows the instantaneous conductance for several different 
temperatures with asymmetry factor of 0.9, fixed Ttot, ^=Q^tot and rectangular DOS. 
Dark blue(dot dot dashed), red(dot dashed), green(dashed), light blue(dotted) and 
black(sohd) curves represent the conductance at T=0.0186, 0.0093, 0.0046, 0.0023, and 
O.OOOQFjot respectively. Red(dot dashed), black(dotted) and dark blue(dashed) curves 
in panel b display the instantaneous conductance vs. time in the long timescale for 
rectangular DOS when the source-drain bias is equal to V=Tk, V=5Tk and V=10Tr- 
respectively with fixed Ttot and 'D=9Ttot for asymmetry factor of 0.9 at T=0.0093rtot- 



4. Discussion 

Our numerical calculations clearly show that the timescale for the decay of the 
conductance oscillations is much larger than the fast timescales set by the couplings 
of the dot level to the leads, Tl or Tl. The timescale does not depend on the width or 
shape of the DOS feature in the leads. It appears that the timescale is related to the 
Kondo resonance. When the Kondo resonance is fully formed at times around 1/Tk, 
the oscillations disappears. 

Further support for the role of Kondo physics in these conductance oscillations 
comes from the observation of a saturation of their amplitude for temperatures and 
source-drain bias below the Kondo temperature (Fig. 4). The effect of temperature 
and source-drain bias on the Kondo resonance in an SET was recently investigated 
in detail [23]. For values smaller than the Kondo temperature, the Kondo resonance 
is fully formed just above the Fermi Level. Increasing the temperature above Tk 
broadens and reduces the magnitude of the resonance. Increasing the source-drain bias 
results in the formation of a split Kondo resonance with strongly reduced intensities. 
For temperatures or source-drain bias well above the Kondo temperature, the Kondo 
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resonance is completely suppressed. Since the frequency of the oscillations is determined 
by the energy difference between the Fermi level and the DOS feature, we believe 
that the oscillations reflect an interference process between the conduction electrons 
associated with the Kondo resonance and the conduction electrons associated with the 
DOS feature. In order to substantiate this explanation we used a simple analytical 
noninteracting Anderson model to calculate the instantaneous conductance following 
the sudden switching of a dot level of a width Fe// equal to that of the final Kondo 
resonance in the interacting model. This effective resonance was switched to a position 
just above the Fermi level of the leads. These calculations were performed numerically 
for both the parabolic and rectangular DOS functions shown in Fig. la and lb. The 
instantaneous conductance was found to display the same conductance oscillations with 
frequencies determined by the bandwidth D as in the original interacting model. The 
oscillations are found to decay over a long timescale determined by Fe//. 

While the Kondo temperature depends only on the total coupling to the leads F^t, 
the total current depends on the asymmetry factor. The reason why the conductance 
oscillations only show up for asymmetric couplings is that the total current through 
the dot is the difference of the currents from the left and right leads Eq. ([6]). The 
interference effect shows up both in the the right and left currents. For symmetric 
coupling, the current oscillations from the right and left leads are out of phase resulting 
in a cancellation of the oscillations in the total current. For asymmetric coupling the 
left and right current oscillations do not cancel resulting in the observed conductance 
oscillations. This result has also been verified numerically for the noninteracting model 
described in the preceding paragraph. 

For dots coupled to leads with very weak DOS features, the conductance oscillations 
will be very small. In Fig. 5 we show the time-derivative of the instantaneous current 
for leads with a parabolic DOS function (Fig. lb) for three different values of the 
bandwidth D. The conductance oscillations here results from the interference of the 
Kondo resonance in the final state and the weak discontinuity in the DOS at the lower 
band edge of the leads. The figure clearly demonstrates oscillations of the instantaneous 
currents with a frequency equal to D in the long timescale. 

The experimental study of the conductance oscillations could possibly be made 
using previously suggested techniques, i.e., by measuring the total charge transport as 
function of pulse duration [15]. For a bandwidth of D=l eV, the oscillation period will 
be of an order of 10"^"^ seconds. Since the current oscillates, electromagnetic radiation 
will be emitted at a rate proportional to [^^]^- For a suitably designed system it may 
be possible to detect the emitted light. For D=l eV the emission would occur in the 
infrared at a photon energy of 0.1 eV. 

5. Conclusion 

In this paper, we employed the time dependent non-crossing approximation to analyze 
the transient current in a single electron transistor in the Kondo regime asymmetrically 
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Figure 5. Time-derivative of the instantaneous current in the long timescale for 
parabolic DOS (Fig. lb). Black(solid), red(dashed) and dark blue(dot dashed) curves 
correspond to D—6.75Ttot, ^—Q^tot and D=13.5rfot for asymmetry factor of 1.0 at 
T=0.0093rtot. 

coupled to two metallic leads with features in their DOS. We show that for asymmetric 
coupling, the conductance can exhibit oscillations which persist for times much longer 
than the timescale for charge relaxation. The origin of these oscillations is an interference 
between the conduction electrons associated with the Kondo resonance and those 
associated with the DOS feature. The amplitude of the oscillations are found to 
depend strongly on temperature and source-drain bias when those exceed the Kondo 
temperature. We hope that our predictions will motivate further theoretical and 
experimental studies. 
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